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熔 盐 堆 流 体 动力 学 模型 降 阶 方法 适用 性 分 析 
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摘要 ”对 于 熔 盐 堆 全 堆 高 保 真 流体 动力 学 计算 ， 即 使 借助 超级 计算 机 的 并 行 计算 能 力 在 面 对 快 速 甚至 
实时 求解 的 问题 仍然 面临 效率 的 巨大 挑战 ， 引 入 和 采用 模型 降 阶 方法 (Reduced Order Modeling, ROM) , 
将 能 够 有 效 解决 这 类 问题 。 基 于 本 征 正 交 分 解 方法 (Proper Orthogonal Decomposition, POD) 与 Galerkin 投 
影 法 ， 引 入 基于 有 限 体积 的 模型 降 阶 方法 (ROM based on Finite Volume approximation, FV-ROM) 和 上 确 
界 稳 定 模型 降 阶 方法 (ROM with supremizer stabilization, SUP-ROMO ， 针 对 液态 燃料 熔 盐 堆 (Liquid Fuel 
Molten Salt Reactor，LFMSR) 层 流 和 满 流 瞬 态 工 况 开展 适用 性 分 析 。 结 果 表 明 ，FV-ROM 方法 在 速度 误差 
和 计算 效率 方面 占有 明显 优势 ， 层 流 和 汕 流 瞬 态 速度 平均 L2 相对 误差 低 于 0.5% 和 0.6%， 且 单 步 长 的 加 速 
比分 别 为 1500 和 1000 倍 左右 ; 相 比 之 下 ，SUP-ROM 方法 在 压力 预测 方面 表现 出 显著 的 优势 ， 层 流 和 满 
流 瞬 态 压力 平均 L2 相对 误差 低 至 0.20% 和 0.38% 。 因 此 ， 通 过 FV-ROM 和 SUP-ROM 两 种 方法 相 结合 的 
方式 进行 熔 盐 堆 流 体 动力 学 速度 场 和 压力 场 预测 ， 能 够 更 加 有 效 地 提高 流体 动力 学 仿真 的 效率 和 确保 瞬 态 
模拟 过 程 计算 可 靠 性 和 精确 度 。 
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Abstract [Background]: For high-fidelity simulations of fluid dynamics in molten salt reactor, even though a 
supercomputer is able to suppress the period of each simulation, the consequent expense is still prohibitively costly. A 
possible way to overcome this limitation is the use of Reduced Order Modelling (ROM) techniques. [Purpose]: Evaluating 
the accuracy of the ROM methods for reconstructing the velocity and pressure fields. [Methods]: Two ROM methods based 
on the Proper Orthogonal Decomposition (POD) with both Galerkin projection, namely FV-ROM (ROM based on Finite 
Volume approximation ) and SUP-ROM (ROM with supremizer stabilization ), are established for fluid dynamics of molten 
salt reactor. Then, both methods are tested on the unsteady cases of liquid-fueled molten salt reactor (LFMSR). [Results]: 
The FV-ROM demonstrates notable advantages in both velocity prediction and computational efficiency. For laminar and 
turbulent transient simulations, the average velocity L2 relative errors are less than 0.596 and 0.696, respectively, with 
acceleration ratios of approximately 1500 and 1000 times for single time steps. Conversely, the SUP-ROM scheme 


demonstrates significant prowess in pressure prediction, achieving remarkably low pressure average L2 relative errors of 
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0.20% and 0.38% for laminar and turbulent transient scenario, respectively. [Conclusions]: The integration of the SUP-ROM 
and FV-ROM for fluid dynamics computations of molten salt reactor could significantly enhances computational efficiency 


and ensure reliability and accuracy of transient simulation. 


Key words Reduced Order Modeling, Molten Salt Reactor, FV-ROM, SUP-ROM, fluid dynamics 
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数值 方法 PC3) ， 即 使 借助 超级 计算 机 的 并 行 计算 能 力 在 面 对 需 要 多 次 求解 ， 快 速 甚至 实时 求解 的 问题 仍然 
面临 效率 的 挑战 。 其 求解 效率 必须 要 有 2-3 个 数量 级 上 的 显著 提升 才能 满足 实际 的 需求 。 模型 降 阶 CROMO 
是 一 种 能 够 保持 精度 的 前 提 下 巨 幅 提升 大 型 工程 计算 效率 的 理论 和 方法 各， 因此 将 ROM 应 用 于 熔 盐 堆 工 
程 设 计 与 分 析 计 算 〈 不 确定 性 量化 、 设 计 优 化 和 运行 控制 ) 具有 重要 意义 。 

模型 降 阶 在 航空 工程 、 造 船 工程 、 汽 车 工程 、 土 木工 程 等 流体 力学 领域 中 应 用 特别 广泛 [0， 在 核能 
系统 领域 近年 来 也 有 相应 的 应 用 。 核 能 系统 的 介入 式 降 阶 方案 中 POD-Galerkin 方法 最 具有 代表 性 ， 如 
German, Alberto 和 Xiang 等 人 利用 POD-Galerkin 方法 解决 带 或 不 带 温度 反馈 的 中 子 扩散 本 征 值 的 计算 问 
AUI: Min 等 人 利用 POD-Galerkin 对 欧洲 压 水 反应 堆 CEPR) 燃料 棒 流 致 振动 特性 进行 分 析 ，Sartori 等 
人 利用 POD-Galerkin 对 铅 快 堆 建立 单 通道 的 耦合 模型 ， 计 算 速 度 提 高 了 至 少 三 个 数量 级 05]，Lorenzi 等 人 
采用 POD-Galerkin 方法 建立 了 铅 冷 快 堆 ALFRED 降 阶 模型 ， 并 模拟 了 ALFRED 的 瞬 态 工 况 09; 陶 阳 等 人 
利用 POD-Galerkin 方法 模拟 铅 铅 快 堆 事 故 过 程 中 上 静 压 室 的 温度 分 布 , 为 确定 热 分 层 机 制 提供 了 一 个 分 析 
T RUT 


SPS ARSE CENE f Fe RAS HE) 这 类 复杂 的 非 线性 系统 ， 对 其 建立 流 场 动力 学 降 阶 模型 还 存在 
两 个 主要 问题 。 首 先 ，Galerkin 投影 到 的 缩减 基 空 间 并 不 能 保留 原 有 的 全 阶 动态 特性 08;， 其 次 ， 由 于 熔 盐 
被 近似 成 不 可 压缩 流体 ， 经 过 投影 后 的 质量 方程 也 具备 了 无 散 度 的 特性 ， 造 成 单方 程 求解 多 变量 的 问题 。 
为 解决 上 述 问 题 ， Lorenzi 等 人 09 提 出 在 降 阶 过 程 中 只 求解 动量 方程 , 忽略 压力 梯度 贡献 的 基于 有 限 体积 的 
降 阶 模型 方法 CFV-ROMD , Ballarin 等 人 29 提出 添加 额外 速度 模 态 来 求解 双方 程 的 上 确 界 稳定 降 阶 模型 方 
ik (SUP-ROMD 。 

为 分 析 和 探讨 FV-ROM 5 SUP-ROM 两 种 模型 降 阶 方法 在 熔 盐 堆 流 体 动 力学 计算 中 的 适用 性 问题 ,本 
文 以 开源 CFD 软件 OpenFOAM 为 基本 框架 建立 熔 盐 堆 流 体 动力 学 全 阶 .FV-ROM 和 SUP-ROM 降 阶 模型 ， 
并 通过 LFMSR (Liquid Fuel Molten Salt ReactonD20 层 流 和 汕 流 瞬 态 工 况 开展 适用 性 分 析 和 评估 。 


1 ”液态 燃料 熔 盐 堆 全 阶 模 型 FOM (Full Order Model) 


由 于 液态 熔 盐 堆 熔 盐 可 认为 不 可 压缩 流体 ， 和 泵 和 换 热 器 处 理 成 多 孔 介质 模型 ， 由 此 可 得 到 质量 及 动量 
方程 : 


V-you=0 (1) 
ð T 
TOL EN (ypu Du)=V-(y(1+n,)| Vu + (Vu) ) 


Wp XP KEV PE, 
其 中 ，p, u, p, y, n, n, DART RE, WE, RA, ÆURK, ONE AE, LEX 


可 观察 到 当 y=1 时 ， 即 为 堆 内 部 流体 区 域 。 忆 ,,, 为 泵 体积 力 ， 为 ( 泵 或 换 热 器 ) 流 阻 ，g, 为 重力 矢量 。 
根据 达 西 速度 w 的 定义 可 表示 为 


uy.) OMA) 6) 
将 式 (3) PAR COD 和 (2) 可 得 到 : 
“eo (4) 
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Pes V (pu ®uy)=V-((n+n, [Vv as 1 


IpsexF og + YF + YPE, 


炊 盐 堆 属 于 闭环 系统 ， 在 边界 上 速度 采用 齐 次 狄 利克 雷 边 界 条 件 Cu, 200 ， 压 力 采 用 诺 伊 曼 边界 条 
件 Vp:n=0。 


(5) 


2 ”POD-Galerkin 方法 
2.1 本 征 正 交 分 解 (POD) 


建立 低 阶 模型 需要 一 组 能 够 在 数值 特征 上 面 代表 原 数 据 系统 的 基 函 数 ， 流 场 内 的 速度 可 以 分 解 成 时 间 
相关 系数 a(t) 乘 以 标准 正 交 空间 基 o" (x) 的 线性 组 合 : 


ua) = Ja Og O) (6) 


获取 的 缩减 基 可 通过 POD (Proper Orthogonal Decomposition Proper) , PGD (Generalized Decomposition 
Proper) , #48 RB (Reduced Basis method) 55:7; 102741, KLH POD 是 一 种 将 一 组 数值 (在 时 间或 
参数 空间 中 ) 压缩 成 一 组 较 少 模 态 方法 ， 这 些 基 能 够 能 捕捉 流 场 中 最 重要 的 信息 。 对 于 稳 态 的 快照 ， 根 据 
Wr xg WIR Eu = (usus... uy) 获取， 快照 总 数 为 Y。 对 于 瞬 态 快照， 需 选 取 个 时 刻 瞬 态 流 场 数据 ， 
因此 其 快照 总 数 为 N, = Nx K ， 可 表示 为 : 


RS uto) eo ees] (7) 


AP, R, 为 速度 样本 。 速 度 的 POD 空间 是 通过 解决 以 下 优化 问题 构建 的 : 


2 


N, 


ile : 
Vpop = argmin =~ > u, 一 S (u, 从 Jew? (8) 
s n=l i=l Po) 
AP, u EHR ROAINTRAR, AEH, RER (8) 等 价 于 求解 如 下 特征 值 问题 : 
CV" S V*A*, (Cy = (usu). (9) 


AP C Fe EE iy PR EB) EA, VERENE A 为 对 角 和 矩阵 ， 其 对 角 元 素 是 特征 值 。 并 通过 
以 下 方式 组 合 形成 POD 基 模 态 ， 并 类 推出 压力 和 消 流 粘度 POD ASP: 


1 x V 
us cbe ; NA 10 
9; Tan Vis 9; hen Peg are ae : (10) 


上 上 式 中 符号 上 标 w,p,n, 代表 了 速度 ， 压 力 和 消 流 粘度 。 根 据 求解 的 特征 值 截断 POD 模 态 来 近似 速度 ， 


uj Gf) s Valuer, PENDO, nD A on)" (11) 


AP a, bue, DANTE, RA Amit GE SEES RB, NNN, aia FN, o 


p? s 


2.2 Galerkin 投影 


将 Navier-Stokes 方程 在 式 (10) 所 得 的 POD 张 成 模 态 空间 上 Galerkin 投影 ， 可 以 将 高 维 的 动力 系统 
转换 成 低 维 的 流动 模型 P729]， 投 影 后 的 式 (5) 下 可 以 写成 以 下 公式 : 
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(V: Pup V 97) 2o) -0 


ð 1 (12) 
| $7 (on 9u,)-V (em [Vun (0) ) Np ors - YF, 7108. e -0 
P(O) 
将 式 C11) WA (12) 中 得 到 : 
Ma=0 
(13) 


— cum N. 
。 T^ 了 万 
pHa+ pa Ca-e Ta-nDa+Pb+ Rls- E ios S 


mw. )- YOR, = 0 


其 中 的 运算 符 如 下 所 示 : 


io alee at e 


Q 


eg) 


pa | (Suy = (g^.8. (j^ pm A 7 (g^ ,yg-r), 
[e] 


ppump,z 


(14) 
式 中 (,) 表示 为 内 积 ，6.() 多 孔 介 质 区 为 1， 非 多 孔 介质 区 为 0。 需要 指出 的 是 ， 作 为 无 散 快 照 的 线性 


组 合 ， 速 度 基 函数 也 是 无 散 的 ， 这 意味 着 矩阵 M 为 0， 导致 无 法 同时 确定 压力 和 速度 模 态 系数 ， 不 满足 
Ladyzhenskaya-Babuska-BrezziB0 条 件 。 为 解决 该 问题 ，Lorenzi 等 人 引入 FV-ROM Z7 TE BI Hs 7] fll iria ifi E 
场 都 遵循 速度 场 的 动力 学 原则 M9， 因此 速度 ， 压 力 ， 满 流 烙 度 场 表 示 如 下 : 


ut) s Ya (tu) Q(x), pGot) & >。 ((,u)g"(Qx), nG.t) s > q(t,u) G(x) (15) 
其 中 基 模 态 可 表示 为 : 
me 1 & 1 & 
Qi? = LuV o == PM, OF = 2g; (16) 
Apo Ai jl Ay j= 


HERRAR (13) 中 ， 只 需求 解 模 态 系数 a (tu) 即 可 。 
Ballarin[?0 等 人 提出 添加 额外 模 态 到 速度 基 模 态 中 的 SUP-ROM 方法 ,该 方法 使 矩阵 M 为 不 为 0， 由 此 
可 获得 两 组 方程 求解 压力 及 速度 的 模 态 系数 。 额 外 模 态 的 快照 是 通过 引入 泊 松 方程 进行 计算 的 : 
As, = -Vp, (= (17) 


其 中 ， 压 力 快照 是 已 知 值 ， 加 入 边界 条 件 〈 5, = 0 ) 即 可 以 求解 。 保 存 上 式 的 快照 R = [sss]: 
根据 POD 方法 求解 特征 值 ， 通 过 式 (11) 计算 并 提取 特征 值 较 大 的 几 个 模 态 ， 将 其 加 入 到 式 (10) 的 速度 
模 态 we 中 ， 使 缩减 基 空 间 中 速度 不 具备 无 散 度 的 特性 。 文 章 中 流 阻 模 态 系数 < 和 汕 流 模 态 系数 e 分 别 通 
过 离散 经 验 插值 和 径 向 基 函 数 方法 获取 的 Bl321。 


2.3 在线 求解 及 误差 统计 
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模型 降 阶 计算 流程 如 图 1 所 示 , 分 为 离线 阶段 和 在 线 阶段 。 离 线 阶 段 包 括 快 照 生 成 和 模 态 获取 等 步骤， 
而 在 线 阶 段 则 首先 执行 降 阶 方程 的 算 子 矩阵 组 装 ， 随 后 求解 代数 方程 组 式 〈13 ) ， 最 后 利用 式 〈11) EW 
流 场 。 


Offline phase Online phase 
Begin 
Samples calculate 
using FOM Assembling ROM 
FV-ROM or SUP-ROM > : solving ROM 
à Reconstruct flow field 
Modes calculate . 
Building reduced : End 
operators : 


DS 


图 1 模型 降 阶 计算 流程 
Fig. 1 Flow chart of construction of ROM 


为 衡量 ROM 的 准确 性 ， 采 用 相对 全 局 L2 误差 作为 指标 ， 因 此 速度 和 压力 误差 可 表示 为 : 


T 


ltron —Urom |Pron — Prom m 


e, MM gs: MR. (18) 


u p 
[trou Po [Prou m 


为 统计 时 间 步 长 上 的 速度 和 压力 相对 全 局 L2 误差 平均 值 ， 定 义 : 


- 1 
e, 二 e, = Dyn (19) 


3  LFMSR 流体 动力 学 降 阶 模型 分 析 


LFMSR 是 一 种 液态 燃料 熔 盐 快 堆 P0， 设 计 的 热 功率 为 800MW， 有 八 个 回路 ， 每 个 回路 由 三 个 主要 部 
件 组 成 : 堆 蕊 、 泵 和 换 热 器 。 一 回路 几何 和 网 格 结构 如 图 2 Br. ded 列 出 了 使 用 的 熔 盐 物性 参数 ， 换 热 
器 孔隙 率 设 定 为 0.5。 泵 作为 动量 源 位 置 放 在 回路 下 方 ， 液 态 人 燃料 流动 的 方向 是 自 下 而 上 ， 网 格 共 划分 了 
44770 单元 。 全 阶 模型 的 计算 是 在 处 理 器 为 Intel@CorerM i5-12500H, 主 频 为 2.50 GHz 的 PC 机 上 进行 ， 本 
文 各 全 阶 算 例 均 采 用 8 核 并 行 计算 。 


or 


q 
E Ur C. 
i ^7 e 
x< 
/ ii [E 
g| w | = 
^| o 
o a 
= || 
人 
wg 
图 2 LFMSR 一 回路 (a) 几何 模型 ， (b) 二 维 网 格 
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Fig.2 LFMSR primary loop (a) geometry, (b) 2D spatial mesh 


表 1 LFMSR 热 物 性 参数 
Table 1 LFMSR Thermophysical 


参数 参数 值 参数 参数 值 
Parameter name Parameter value Parameter name Parameter value 
密度 热膨胀 系 关 
idu 3 4200 TERM -— 4.60E-4 
Density / kg:m* Thermal expansion coefficient/K 
热 容 3 Bap a 

ind Por, 950 Binet 2.945 
Specific heat/ J:kg +K Prandtl number 


Pump Power (%) 
Pump Power(%) 


T T - T » T = T T 
0 50 100 150 200 250 300 350 0 1 2 3 4 5 
Time (s) Time(s) 


图 3 瞬 态 工 况 下 功率 变化 (a) 层 流 ， (b) mi 


Fig.3 Relative pumping power of the core for transient problems (a) Laminar, (b) Turbulence 


图 3a 展示 了 0-350s 层 流 工 况 下 泵 功率 变化 曲线 ， 初 始 时 刻 泵 动量 源 设 置 成 10Pa/m， 每 阳 5 秒 获 取 一 
个 快照 共 获 取 TO 个 快照 。 图 3b 展示 了 0-5s 灌流 工 况 下 泵 功率 变化 曲线 ， 初 始 时 刻 泵 动量 源 设置 成 
20000Pa/m， 每 隔 0.05 秒 获 取 一 个 快照 共 获取 100 个 快照 。 图 4 为 层 流 和 湛 流 工 况 下 压力 、 速 度 、 流 阻 和 
SUP 的 模 态 特征 值 曲线 ， 可 以 观察 到 层 流 工 况 下 各 流 场 模 态 能 量 衰减 较为 迅速 ， 并 且 前 10 模 态 基本 包含 
流 场 的 大 部 分 信息 。 瀚 流 工 况 下 满 流 粘度 和 速度 衰减 较 慢 ， 因 此 需要 选择 更 多 的 模 态 才能 反映 流 场 信息 。 
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图 4 速度 ， 压 力 和 流 阻 的 归 一 化 模 态 特征 值 曲 线 (a) EX. (b) du 


Fig.4 The normalized eigenvalues of modal for velocity, pressure and flow resistance (a) Laminar, (b) Turbulence 
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FV-ROM fi 


对 比 了 5 
都 呈现 了 一 利 


Velocity relative L2error(%) 


Pressure relative L2error(%) 


200 300 
Time(s) 


100 200 300 
Time(s) 


图 5 层 流 工 况 下 不 同 模 态 下 FV-ROM 和 SUP-ROM If] L2 相对 误差 (a) 速度 ， (b) 压力 
Fig.5 L2 norm relative error of FV-ROM and SUP-ROM in different modes at laminar scenario (a) Velocity, 


降 阶 模型 的 计算 量 通常 与 所 采 


JA BLAS CR 


Il SUP-ROM 方法 在 使 用 不 同 模 态 数 下 的 性 能 表现 显 
模 态 数量 下 的 L2 相对 误差 ， 


4 


F 8 时 速度 误差 基本 重合 ， 


不 稳定 的 趋势 。SUP-ROM 对 比 了 5 组 模 态 数 ， 除 了 第 一 组 所 选 模 态 累积 能 量 占 比 过 低 导 致 预 


得 尤为 重要 。 图 5 展示 了 如 
中 Nu. Nsup 和 Np 分 别 表示 所 选择 的 速度 .SUP AU Fs 77 RAS CR 
速度 模 态 数 。 当 速度 模 态 数 大 了 


(b) Pressure 


成 正比 ， 但 不 意味 着 计算 精度 会 随 之 提升 ， 基 


此 ， 评 估 


E 层 流 工 况 下 不 同 


E 


Eo 


速度 模 态 数 越 大 压力 误差 有 所 减少 ， 但 


FV-ROM 


测 的 速度 误差 较 大 ,后 四 组 随 着 SUP 模 态 数 增 大 速度 的 平均 误差 和 最 大 误差 也 在 扩大 ,其 主要 原因 是 添加 


了 额外 SUP 模 态 而 偏离 了 原 有 的 物理 
Sb 发 现 ， 在 使 用 最 大 模 态 组 时 ， 虽 然 在 200s 附近 压力 会 H 


REUS, SUP 模 态 的 增多 反而 对 速度 的 预测 造成 不 条 


误差 和 最 大 误差 的 角度 来 看 ， 这 三 组 模 态 数 之 间 的 差异 并 不 显著 。 


| 的 影响 。 通 过 图 
时 现 一 个 非常 小 的 误差 ， 仅 为 0.005%， 但 从 平均 


表 3 为 不 同 模 态 下 速度 和 压力 平均 L2 相对 误差 及 最 大 L2 相对 误差 ， 通 过 综合 对 比 ， 本 文选 取 的 


FV-ROM 的 Nu 为 10，SUP-ROM 为 Nu=10、Nsup=8 和 Np=6， 并 将 两 者 的 的 结果 进行 对 比分 析 。 对 于 速 
度 场 预测 ，SUP-ROM 平均 误差 为 1.037% 高 于 FV-ROM 的 0.444%， 最 大 误差 分 别 为 3.704% 和 1.110%。 对 


ES 


表 3 速度 和 压力 平均 L2 相对 误差 = 及 最 大 L2 相对 误差 e.,、 


Table 3 Velocity and Pressure average L2 norm error and max L2 norm error 


压力 场 的 预测 ， 从 图 Sb 中 可 以 直观 看 到 SUP-ROM 的 显著 优势 ， 误 差 保持 在 0.78% 以 下 ，FVROM 只 有 
些 时 刻 点 上 预测 精度 接近 于 SUP-ROM。 


Wr POD 模 态 数 

ee d POD modes El Gy max (%) e,l €, max (%) 
Nu-4 0.904/2.021 5.932/22.362 
Nu=6 0.504/1.164 2.941/14.65 

FV-ROM Nu=8 0.448/1.117 1.482/12.176 
Nu=10 0.444/1.110 0.865/10.046 
Nu=12 0.444/1.124 0.747/8.985 
Nu=4 Nsup=2 Np=2 1.385/3.985 0.241/0.579 
Nu=6 Nsup=4 Np=4 0.982/3.206 0.197/0.753 

SUP-ROM Nu=8 Nsup=6 Np=6 1.034/3.698 0.197/0.745 
Nu=10 Nsup=8 Np=8 1.037/3.704 0.187/0.778 
Nu=12 Nsup=10 Np=10 1.048/3.778 0.188/0.778 


前 文中 分 析 的 L2 相对 误差 只 是 全 局 指标 ， 我 们 给 


XXXXXX-7 


了 速度 和 压力 在 瞬 态 模拟 结束 350s 时 的 绝对 误差 


Ei dX 术 2OXX, XX: XXXXXX 


如 图 6 所 示 。 可 以 看 出 , 在 350s It, FV-ROM 和 SUP-ROM 预测 的 最 大 的 速度 误差 都 出 现在 堆 芯 出 口 处 附 
近 ， 两 者 的 速度 误差 分 布 基本 类 似 ， 但 FV-ROM 的 最 大 绝对 误差 远 远 小 于 SUP-ROM. FV-ROM 较 大 的 压 
力 误 差 集中 在 泵 和 换 热 器 区 域 附 近 ， 而 SUP-ROM 在 堆 芯 出 口 处 且 范 围 小 。 
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6 350s 时 FV-ROM 5 SUP-ROM 速度 和 压力 绝对 误差 


(a) 速度 (FV-ROMO , (b) 速度 CSUP-ROMD (c) 压力 CFV-ROMO , (b) 压力 CSUP-ROMD 
Fig.6 Absolute error of velocity and pressure field between FV-ROM and SUP-ROM at 200s 
(a) Velocity (FV-ROM), (b) Velocity (SUP-ROM), (c) Pressure (FV-ROM), (d) Pressure (SUP-ROM) 


[+| * FV-ROM 
* SUP-ROM 


10000 F 


1000 Ed 


Acceleration Ratio(FOM/ROM) 


100 H 1 1 1 1 i 
0 50 100 150 200 250 300 350 


Time(s) 


图 7 单个 步 长 时 间 加 速 比 
Fig.7 Single step time acceleration ratio 
图 7 为 整个 层 流 瞬 态 工 况 下 每 个 时 间 步 FV-ROM 和 SUP-ROM 的 加 速 比 。FV-ROM 比较 稳定 在 1500 
倍 左右 ， 总 体 上 其 加 速 比 是 优 于 SUP-ROM。SUP-ROM 波动 较 大 ， 主 要 原因 是 SUP-ROM 在 泵 功率 平稳 阶 
BIE A SAGE BE BEER 


32 WALAA ROAT 
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Velocity relative L2error(%) 


—sUP-ROM (Nu-5 Nsup-2 Np=2} 


1000 E 


E|—— FV-ROM (Nu-20) 
[—— Fv-RoM (Nu=25) 


) 


Pressure relative L2error(%) 


图 8 


T T T T 
if 2 3 4 


Time(s) 


Time(s) 


MALI BAS BLAS F FV-ROM 和 SUP-ROM 的 L2 相对 误差 (a) 速度 ， (b) 压力 
Fig.8 L2 norm relative error of FV-ROM and SUP-ROM in different modes at turbulent scenario (a) Velocity, (b) Pressure 


为 了 评估 在 满 流 工 况 下 FV-ROM 与 SUP-ROM 方法 使 用 不 同 模 态 的 表现 ， 同 样 对 比 5 组 模 态 数 ， 其 速 
度 和 压力 的 L2 相对 误差 如 图 8 所 示 。 观 察 到 FV-ROM 速度 模 态 数 为 15 时 速度 误差 基本 保持 在 1% 以 下 ， 


| 


但 都 低 于 1.0396. 


当 模 态 数 大 于 20 时 速度 误差 有 所 增 大 ， 压 力 在 3.5s 后 出 现 明显 的 偏差 ， 
SUP-ROM 采用 模 态 


组 合 为 (Nu=15 Nsup=3 


误差 最 高 达到 100% 以 上 。 当 


Np=3) 时 ， 速 度 预测 是 最 准确 的 ， 压 力 误差 略 高 于 其 他 模 态 组 


0.387%， 最 大 的 误差 不 超过 1.1%，SUP-ROM 表现 出 显著 的 优势 。 
4 ”速度 和 压力 平均 L2 相对 误差 z 及 最 大 L2 相对 误差 e,、 


Table 4 Velocity and Pressure average L2 norm error and max L2 norm error 


表 4 为 不 同 模 态 下 平均 L2 相对 误差 及 最 大 L2 相对 误差 ， 通 过 综合 对 比 ， 本 文选 取 的 FV-ROM 的 Nu 
为 15，SUP-ROM 为 Nu=15 Nsup=3 Np=3， 并 将 两 者 的 的 结果 进 和 
平均 L2 相对 误差 为 0.599% 低 于 SUP-ROM 的 0.914%。 在 压力 的 预测 上 ，SUP-ROM 平均 L2 相对 误差 为 


于 对 比分 析 。 在 速度 预测 上 ，FV-ROM 的 


降 阶 方法 POD 模 态 数 zd - , 
ROM method POD modes €,! €, max (%) e, £p max ( o) 
Nu=5 1.008/1.993 27.595/43.227 
Nu=10 1.085/2.312 36.628/58.602 
FV-ROM Nu=15 0.599/0.989 13.217/27.148 
Nu=20 0.933/1.939 30.938/155.336 
Nu=25 0.920/1.929 29.935/50.816 
Nu=5 Nsup=2 Np=2 1.075/1.847 0.406/1.071 
Nu=10 Nsup=2 Np=2 1.007/1.679 0.263/0.945 
SUP-ROM Nu=15 Nsup=3 Np=3 0.914/1.510 0.387/1.027 
Nu=20 Nsup=4 Np=4 0.959/1.636 0.259/0.938 
Nu=25 Nsup=5 Np=5 0.934/1.576 0.253/0.933 
速度 和 压力 在 瞬 态 模拟 结束 5s 时 的 绝对 误差 如 图 9 所 示 。 可 以 看 出 , FV-ROM 和 SUP-ROM 预测 的 较 


大 的 速度 误差 都 集 


的 压力 误差 都 集中 在 和 泵 和 换 热 器 区 域 ，SUP- 


在 堆 蕊 下 方 ， 并 且 范 围 


于 压力 的 预测 是 比较 准确 的 。 


不 广泛 ， 两 者 的 速度 误差 分 布 基 
ROM 则 是 在 堆 芯 上 部 的 极 小 的 
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9 5s 时 FV-ROM 5 SUP-ROM 速度 和 压力 绝对 误差 


(a) 速度 CFV-ROM) , (b) 速度 (SUP-ROM) (c) 压力 (FV-ROM) ， (b) 压力 (SUP-ROM) 


Fig.5 Absolute error of velocity and pressure field between FV-ROM and SUP-ROM at 5s 


(a) Velocity (FV-ROM), (b) Velocity (SUP-ROM), (c) Pressure (FV-ROM), (d) Pressure (SUP-ROM) 


* FV-ROM 
* SUP-ROM 


10000 F 


1000 


Speedup rati(FOM/ROM) 


Time(s) 


10 单个 步 长 时 间 加 速 比 


Fig.10 Single step time acceleration ratio 


10 ARN Tis V BE AS IL EP EIN [2 FV-ROM 和 SUP-ROM 加 速 比 ，FVROM 加 速 比 集中 在 1000 
AA, SUP-ROM 总 体 加 速 比 在 800 左右 低 于 FV-ROM， 其 原因 是 SUP-ROM 组 装 的 算 子 矩阵 较 大 ， 计 算 
相对 耗 时 。 上 述 表明 了 在 湛 流 瞬 态 工 况 下 即使 全 阶 模型 采用 了 八 核 并 行 计算 , 降 阶 后 加 速 比 还 是 比较 理想 。 


4 ”结语 


针对 熔 盐 堆 流体 动力 学 计算 效率 问题 , 引入 了 FV-ROM 和 SUP-ROM 两 种 模型 降 阶 方法 , 并 以 LFMSR 
层 流 和 汕 流 的 瞬 态 工 况 作为 测试 案例 ， 详 细 分 析 了 两 种 方法 的 流 场 和 压力 场 预测 性 能 和 计算 效率 。 数 值 计 
算 结 果 表 明 ， 两 种 方法 都 能 够 显著 提高 计算 效率 ， 其 中 FV-ROM 在 适当 的 模 态 下 能 够 提供 更 加 精确 的 速度 
信息 ， 而 SUP-ROM 则 有 助 于 更 加 准确 地 捕捉 压力 场 的 变化 。 因 此 ， 通 过 结合 这 两 种 方法 ， 能 够 在 保证 燃 
盐 堆 流体 动力 学 计算 可 靠 性 和 精确 度 的 情况 下 ， 显 著 提 高 流体 动力 学 仿真 效率 ， 为 熔 盐 堆 流 体 动力 学 计算 


分 析 提 供 高 效 的 理 * 


企 和 方法 


o 


作者 贡献 声明 林 铭 负责 建立 全 阶 及 降 阶 模型 ， 数 据 分 析 及 文章 撰写 : FR TTR OTTER, OEC 
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